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1 Introduction 

In this paper, we introduce a new representation of many body electron wave function 
and a few calculation results of the ground state energies of many body systems using 
that representation, which is systematically better than the Hartree-Fock approximation. 

The fundamental principle of condensed matter physics and chemistry is given in the 
many body schrodinger equation, which is 

N 

H = E("(1/2)V, 2 + v(r t )) + £ l/|r, - rj \, (1) 

Hip(xi ■ ■ ■ xn) = Eip(xi ■ ■ ■ xn). (2) 

This is a schrodinger equation of N electrons in the external potential v(r) from nuclei 
under the Born-Oppenheimer approximation. The electron wave function if)(xi ■ ■ ■ xn) 
must be antisymmetric in the exchange of arbitrary two spin coordinates Xi and Xj. One 
of the central problem in condensed matter physics and chemistry is to find the solution 
of ([2]) and the ground state energy E of the given system. There are many attempts to 
solve this problem. In this section, we briefly look over the Hartree-Fock approximation 
(HF) and the Density functional theory (DFT) or Local density approximation (LDA). 

In the Hatree-Fock approximation, the many body wavefunction ip(xi ■ ■ ■ xn) is ap- 
proximated by a single slater determinant ipsLi and this gives limitation to the accuracy of 
calculated ground state energy Ehf- The reason is that an antisymmetric wave function 
if)(xi ■ ■ - xn) is not always given as a single slater determinant ipsL- Rather, ip(xi ■ ■ -x^) 
can be expanded in a linear combination of mCn slater determinants in the space of given 
orbital set ipi, • • • ,i()m- in this way, one can rearrange the many body electron prob- 
lem into the diagonalization of mCn x m Cn matrix. This method is known as Full CI 
method. Full CI always gives the exact ground state energy in the space of a given orbital 
set ipi, ■ ■ ■ ,ipM- However, when the number of the electrons N and the orbitals M in- 
creases as much as a few dozens, the dimension of the matrix mCn increases exponentially 
with N and M and Full CI calculation becomes practically impossible. 

In the Hartree-Fock approximation, the variables of the variation are the N orbitals 
ipi, ■ • ■ , tpjf. In contrast, in the density functional theory, the variable is the one electron 
density p(r). This p(r) can be uniquely deduced from the given antisymmetric wave 
function ip{x\ ■ ■ •x^f). The essence of the DFT is that the expectation value of the sum 
of the kinetic energy T and the coulomb repulsion energy U in the ground state can be 
given by a unique functional of the one electron density, E[p(r)]. However, the exact 
form of this functional F[p(r)\ is not known until today. Therefore, in LDA calculation, 
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several approximated form of this functional is used and they are not exact. In LDA, 
the v-representability of p(r) is generally assumed, if one uses integers for the occupation 
numbers. In this case, the calculated ground state energy E LD a and the one electron 
density p(r) are derived from a non interacting single slater determinant ipsL- This may 
be the reason why the LDA calculation does not work well in the so called strongly 
correlated systems. 

In the next section, we introduce a new representation of the antisymmetric wave 
function which is an extension of the slater determinant and therefore not always non 
interacting. 



2 A new representation of many body wave function 

In the many body schrodinger equation (j2J), the variable is the A body antisymmetric wave 
function ip{x\ ■ ■ - xn)- However, this function is apparently not suited for any variational 
calculation because if one takes m spatial grids for one variable Xi, then the total grids of 
the function is proportional to m N and becomes progressively impossible to stock in the 
memory of the computer when he increases A to only a few dozens. One way to express 
this wave function ip(xi ■ ■ ■ xjf) is to expand it in a given M orbital set • • • , ipM- In 
this way, the wave function is expressed as 

ip(xi ■ ■ - x N ) = Ai-iN^hixi) ■ ■ ■ ip iN {x N ). (3) 

ivijsr 

In order to keep ip(xi ■ ■ -x^) antisymmetric, the rank-A and dimension-M tensor A^...^ 
should be antisymmetric tensor. This antisymmetric tensor has mCn degree of freedom 
and still be hard to take variation except for very small systems. 

Here we introduce a new representation of this antisymmetric tensor A^...^: 

K 

Ai v ..i N = ^ C i C ix ' ' ' C ijv e *l*2 e «l«3 ' ' ' e iN-liN- (4) 

1=1 

In this representation, each c},- ■ ■ , cf represents a vector of dimension M. For this reason, 
let us tentatively call this representation as the vector product, ci, • • • , ck are arbitrary 
coefficients. e»j is a rank-2 and dimension-M tensor which is defined by 



1 (i < j) 

e i:j = { -1 (i>j), (5) 
(i=j). 

It is easy to verify that e^...^ = ■ ■ ■ ej JV _ 1 i JV , A(A — l)/2 product of ey, is rank-A 
and dimension-M antisymmetric tensor which takes the value of only 1, —1 or for any 
indices %\ ■ ■ -i^. Here we give the proof. If A = 2, it is obvious from the definition that 
e ili2 is antisymmetric. Let us assume that e^...^ = e ili2 ■ ■ ■ e iN _ liN is now antisymmetric 
and consider about the rank-A + 1 tensor 

We only need to prove that this tensor is antisymmetric for indices z'i and zjv+i- One can 
rewrite the above form to 
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When one exchanges the indices i\ and ijv+i, tensor e^...^ does not change. Tensor 
^hi 2 ' ' ' e hi N an d e «2jjv+i ' ' ' e «jvijv+i change into each other with the same factor (— l)^ -1 , 
so their product does not change. And finally, the tensor e iliN+1 changes its sign by — 1. 
Therefore, the tensor 6j 1 ...i Ar+1 changes its sign by —1 in this operation. Thus the argument 
is proved. It is trivial from the definition that the value of the tensor e^...^ takes only 1, 
-1 or 0. 

From this, one can conclude that the tensor A^...^ is antisymmetric. We have proved 
that e^...^ is antisymmetric. The vector product c\ ■ ■ ■ c\ is apparently symmetric. 
Therefore, c$ • • • c\ e^...^ is antisymmetric. A il ... iN is a linear combination of these anti- 
symmetric tensors for % — 1, • • • , K, then antisymmetric. 

Notice that the representation (j4]) is an approximative form of general antisymmetric 
tensor. Antisymmetric tensor Ai x ...i N has mCn elements, and belongs to the dimension 
mCn vector space. Therefore, if K =m Cn and all c- x ■ ■ • c\ e iv .. iN are linearly independent 
for i = l,---,K, then any Ai 1 ...i N can be expanded in the form of (jlj). However, in the 
following section we will see that one can well approximate at least some antisymmetric 
tensors which are solutions of the many body problem, with the condition K = 1. 

Apparently, the vector product for K = 2 is an extension of the vector product for 
K — 1. Generally, the vector product for K + 1 is an extension of the vector product for 
K. Here we see that the vector product for K = 1 includes all slater determinants for a 
given orbital set ipi, • • • , ipM- When K = 1, the vector product is given as 

A^-.-ij^ CCjj • • • Ci N €{ 1 i 2 ■ ■ ■ ^i N _ x i N - (8) 

We only need to consider the case which ii = 1, ■ ■ ■ = N. When it is possible to expand 
the wave function ijj(xi ■ ■ ■ xn) in the orbital set ■ ■ ■ , ipN, the antisymmetric tensor 
Ai v ,,i N has degree of nCn = 1 and thus uniquely determined under the antisymmetric 
condition, with an arbitrary factor. The slater determinant for ipx, ■ ••,'0iv is of course 
antisymmetric. Therefore, we only need to show that the vector product for K = 1 can 
be antisymmetric when expanded in orbitals ipi, - ■ ■ , ipN and zero when expanded in other 
orbitals. This condition is satisfied when one takes the vector for 

c ={ 1 (*i = V"-,A0, , 9) 
n \ (others). 1 ' 

Thus the argument is proved. 

It is not difficult to see that the vector product with K = 1 (ipy P ) is indeed an extension 
of the slater determinant if the number of the orbital M is larger than the number of 
electrons N. In other words, there are wave functions which are representable in the form 
of tpyp but not in a slater determinant. The characteristic of the slater determinant is 
that the first order reduced density matrix (71) of the slater determinant has the same 
eigenvalues 1/N for iV natural orbitals. One can easily see that 71 of an arbitrary wave 
function which is representable as ip VP has generally different eigenvalues A.; < 1/N by 
simply take random values for the vector and calculate 71 and its eigenvalues. This 
fact indicates that ipyp is generally interacting and not always representable by a single 
slater determinant. 

As we discussed above, the vector product with K = 1 is an extension of a single 
slater determinant. For this reason, one can obtain the same or lower ground state energy 
Eq by using as variational parameters, compared to the Hartree-Fock approximation. 
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This is possible when one uses a orbital set ipi, - ■ ■ ,ipM which includes the Hartree-Fock 
derived orbitals • • • ,ipN- One can obtain further lower energy by the vector product 
method when he takes variation also for orbitals ipx, ■ ■ ■ , ipM- This variation of orbitals is 
possible in various way. One way is to take unitary transform of these orbitals and take 
variation for the elements of the unitary matrix. The details of this transform are to be 
explained in the following section. 

3 Energy calculation in the vector product method 

Generally, the expectation value of the energy E of the N electron system in the normal- 
ized state ip(xi, • • • , xn) is given as 

E — J ■ ■ ■ J dxi ■ ■ ■ dxNip*(xi, ■ ■ ■ , xiy)Hip(xi, ■ ■ ■ , xn). (10) 

However, the hamiltonian of electrons is generally two body and in the following form: 

H = E(-(1/2R 2 + «(r<)) + E Vh - ^1, (11) 
Therefore, one can obtain E from the normalized second order reduced density matrix 

(72) m 

E = N Jdx 1 [(-(l/2)Vj + v(r 1 )) ll (x' 1 ,x 1 )} x[=Xl 

+ N(N-l)/2f fd Xl dx 2 - n 2 (xix 2 ,x 1 x 2 ), (12) 

J J \V\ — v 2 \ 

72 «z< 2 , IlX2 ) J...J.,r.- dx.nx^x, • • • XN ),^ lX2X3 ...«„), (13) 

TlM,*) -/*»»««».«.*»). <") 
One can obtain the tensor representation of 72 by substituting (Ej) for (|T3l) : 

"( 2 {x 1 X 2 ,X Z X 4 ) = E C iii2i 3 iA^t{ x Mi 2 Mfas (^3)^4 (^4), (15) 
= E Ai^i' 2 i 3 ...i N Ai 1 i 2 i 3 ...i N . (16) 

Here, orbitals • • • , ipM are assumed to be orthonormal. One can write down the tensor 
representation of E as following: 

E = hi 1 i 2 i 3 i 4 Ci 1 i 2 i 3 i 4 /n, (17) 

n = (ip(x 1 ---x N )\ip(x 1 ---x N )) 

— /] (18) 

*1«2 

^*1«2*3*4 = ^\Pil%3 ^1113)^2*4 7> ^2112*3*45 (19) 



t hi3 = I dx 1 ^ 1 {x 1 ){-l/2)V 2 Mx 1 ), (20) 
dxi^j* {xxjv^i)^ (xx), (21) 

w ili2 i 3 i 4 = J J dx l dx2^* 1 {xx)ip* 2 {x 2 ) , _ ^ M )i 3 {x 1 )ipi A {x2). (22) 

Next, we will see the representation of 72 in vector product. Here we assume real 
values for the components of the vector product q and d a for simplicity. We need to 
substitute (HI) for (fT6l). The result is 



K 

5-/ C J3 C j'3 e *l-?'3 e «2j3 e «3i3 e *4j3 ' * ' C j N^j N^hj N^i23 N^hj N^iij N^ji—j N 
3S---3N 
K 

= C « C j C ii C i 2 C i3 C i4 e H«2 e j3M ' ^1^21314,1 (23) 

d- • = e 2 • 

a 3:y-3N — ' j::-J\ 

f 1 (j3 ■ ■ • jiv : all different), . . 

\0 (others). 1 ' 

There is a method to calculate the tensor Illi 2 i 3 i 4 - When one omits the indices ilizizib 
and ij, Illi 2 i 3 i 4 has a following form: 

/ = 22 dj 3 ■ ■ ■ dj N dj 3 ...j N , (25) 

33— 3N 

a 3k = C j k C j k e iljk e i2jk e i3jk e U3k- (26) 

From the definition of the tensor dj 3 ...j N , one can conclude that the sum in (f25|) is taken 
over all permutation (j 3 • • • j^) f or jz • ■ ■ jjv = 1, • • • , M with the condition that all j$ ■ ■ ■ ]n 
are different. Therefore, the value of / in (f25l) can be represented as the M — (N — 2)-th 
order coefficients of the following polynomial: 

f(t) = (N-2)\(t + a 1 )---(t + a M ). (27) 

Then the problem is how to calculate the M — (JV — 2)-th order coefficients of a given 
polynomial f(t). Notice that now f(t) is an M-th order polynomial of t. Then f(t) can 
be expanded in a following form: 

f(t) = b + b 1 t + --- + b M t M , (28) 

I = b M -(N-2). (29) 

One can derive the M — (N — 2)-th order coefficients &M-(iV-2) by solving a linear algebra 
problem. One can define a vector b of M + 1 dimension from the unknown coefficients 
bo, ■ ■ ■ , 6m- One can have arbitrary M + 1 different values for to, • • ■ , t m ■ Then he has 

/(to) = 6 + Mo + --- + &mC 

= &o + Mm + --- + Mm- (30) 
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One can solve this linear algebra problem with given values for £ > ■ ■ ■ , tu and / (to) /(^m) 
and find all the coefficients fej. To do that, one only need to calculate the inverse matrix 
of a matrix = t\ and multiplies it to a vector composed of f(t ), • • • , f(t M ). In this 
way, one can calculate the value of &m-(jv-2)> an d therefore the value of tensor / for each 
indices i\i 2 ij,i^ij. Then one can calculate the matrix element Cj^^ from c$ and d v 

One has to take calculation steps proportional to M 2 for each calculation of Z^^, 
because he has to calculate the value of M-th order polynomial f(t) for about M times. 
To calculate Cj^^, one has to take this M 2 step for M times for each indices, therefore 
he needs calculation steps proportional to M 6 for each calculation of the matrix C^^. 



4 Variation of orbit als 

In this section, we will consider about the variation of orbitals used in the expansion of 
many body wave function. The result of the vector product method is depending of the 
choice of orbitals, then one has to take variation for orbitals in the vector product method 
in order to obtain lower energy. Let us assume that we use orbitals ipi, - ■ ■ , ipu and a basis 
set 0i, ■ ■ • , 4>b- Here, B > M and orbitals are expanded in the basis set: 

Ux) = J2U lJ <j )3 (x). (31) 
i=i 

We assume that the basis set is orthonormal. Then, dimension- B vectors Uij, • • ■ ,Umj 
should be orthonormal in order to maintain orbitals ipi, ■ • ■ , ipM to be orthonormal. Next, 
we will see how this expansion is applied on the expression of the energy matrix hi 1 i 2 i :j i i 
appeared in ffTTl) . The definition of hi li2 i 3ii is 

h hh i 3ii = J J dx 1 dx 2 ip il (x 1 )i/j i2 (x 2 ) [iV(-(l/2)Vi + u(t*i)) 

+ N{N ~ 1 \ 1 (32) 
2 \r 1 -r 2 \ 

By substituting (jHTl) for f[3"2"j) . one gets following expression: 

B 

^ili2i3M — ^iih^iih^^h^iAOi-^hhhH-i (33) 

Jlj2j3j4=l 



J J dx 1 dx 2 (f) jl (x l )(f) j2 (x 2 )[N(-(l/2)Vl + v(r 1 )) 
N(N-l) 1 

+ -i Mj^x^x,). (34) 

/ |t*i — r 2 \ 

This matrix h ili2i3ii can be calculated by taking proportional to MB A steps from the 
fixed matrix Hj x j 2 j 3 j A . The variation of orbitals is possible when one takes variation for 
the matrix Uij with maintaining vectors U\j, ■ • ■ ,UMj to be orthonormal. Notice that 
when B = M, Uij is a unitary matrix or orthogonal matrix. 
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5 Results of calculation 



In this section, we will report on the results of the calculation of ground state energy of 
carbon ( 6 C) and oxygen ( 8 0) atom, using the vector product method with K = 1. The 
author of this paper admits that he is not a specialist of computational physics. He does 
not emphasize the numerical accuracy of the results shown in this section. The purpose 
of this section is to show that the calculation of the expectation energy with the vector 
product is possible and one can obtain better results with the vector product compared 
to the Hatree-Fock. For calculations of the diagonalization in Hatree-Fock and Full CI, 
we used CLAPACK which is LAPACK usable in C. 

First, we will report on the results of 6 C (E V p) by using the vector product method 
with K = 1 with 14 atomic orbitals as a basis set, which is Is 2 ,- ■ -, 3s 2 and 3p 2 , with the 
effective nuclear charge Z = 5.5. In this calculation, we use M = 14 orbitals and set B 
to 14. We took variation for dimension-14 vector c\ and a 14 x 14 orthogonal matrix 
Uij. The expectation value of the energy E is given explicitly as a polynomial of cj and 
Uij. Therefore, we can explicitly define the differential value of E for each cj and U^. 
Then, we took variation for c£ and Uj by using the steepest descent method. We started 
from randomly chosen value for cj and a unit matrix for Uj as an initial value. We also 
calculated the ground state energy in this basis set by the Hartree-Fock method (Ehf) 
and Full CI method (Epci)- The results are following: 



M=14 6 C 


Ehf 


-37.661 


M=14 6 C 


Evp 


-37.688 


M=14 6 C 


Efci 


-37.708 



In this results, one can see that even though we set K = 1 and used only one 4>vp as m any 
body wave function, he can obtain well better results compared to the Hartree-Fock. If 
one takes K larger than 1, then it is sure that he will obtain further better results. This 
result Eyp = —37.688 is obtained by the steepest descent method and not yet converged. 
Therefore we may rather say Evp < —37.688. 

Next, we will report on the results of 6 C and 8 atom using the vector product 
method with K = 1 (Eyp) with 19 atomic orbitals as a basis set, which is Is 2 ,- • •, 3p 6 and 
3dg 2 2_ r 2. We set the effective nuclear charge as Z = 5.5 for 6 C and Z = 7.5 for 8 0. We also 
calculated the ground state energy in this basis set by the Hartree-Fock method (Ehf)- 
We also compare these results with results in a literature E HF , E% FT and E EXP ^\, which 
are results in Hartree-Fock, Density functional theory and experiment, respectively. The 
results are following: 



M=19 6 C 


Ehf 


-37.666 


M=19 6 C 


Evp 


-37.790 


6 c 


^HF 


-37.702 


6 c 


F a 

J^DFT 


-37.479 


6 c 


E a 

^EXP 


-37.858 




M=19 8 


Ehf 


-74.908 


M=19 8 


Evp 


-74.956 


8 o 


Ehf 


-74.858 


8 o 


rpG 

^DFT 


-74.532 


8 o 


Eexp 


-75.113 
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Here, Eyp is not converged for both 6 C and gO. From these results, one can conclude 
that the results of the vector product method with K = 1 can be better than the results 
of Hatree-Fock. If one takes K larger than 1 then he will obtain further lower results for 

Eyp. 

In the calculation of the vector product method, one needs proportional to M 6 steps 
for an evaluation of the matrix 72 . In , There are M variables for variation. Therefore 
we took proportional to M 7 steps for the variation of matrix 7 2 . In the variation of 
the orthogonal matrix Uij, one needs to take proportional to M 5 (MB 4 ) steps for each 
transformation of the energy matrix h ili2 i 3 i 4 . We took variation for each rows of the matrix. 
Therefore we spend proportional to M 6 (M 2 B i ) steps for the variation of orbitals. Then 
the estimated calculation time T for the vector product method is 

T ~ 0(M 7 ) + 0(M 2 B 4 ). (35) 

The total amount of the calculation time for Eyp was about a few hours for qC and sO 
(M = 19) with a ~ 3GHz CPU. In our calculation, we spend larger amount of time for 
the term 0(M 6 ) = 0(M 2 B 4 ) compared to 0(M 7 ). 

6 Conclusion 

We propose a new representation of many body electron wave function, namely the vector 
product. We also propose its application as a post Hatree-Fock method to evaluate the 
ground state energy of many body electron systems. The results of the vector product 
method will converge to the results of the Full CI method when one takes sufficiently large 
value for the parameter K and keeps each vector products ipy P linearly independent. We 
obtained systematically better energy results compared to the results of the Hatree-Fock 
method for qC and sO atoms. The estimated calculation time T for the vector product 
method is T ~ 0(M 7 ) as a function of orbital number M. It is expected that using higher 
spec CPUs, one can obtain the ground state energy better than Hartree-Fock in more big 
systems. In the vector product method, one can simultaneously variate the orbitals used 
in the calculation. This is an advantage of the vector product method compared to other 
methods such as CI method in which orbitals are fixed during the calculation. In the 
vector product method, one can simultaneously obtain the many body wave function of 
the system. This means that one can obtain many physical quantities of the system at the 
ground state. For example, the off diagonal long range order (ODLRO) in solid which is 
related to the superconductivity can be explained in a way that the maximum eigenvalue 
of the second order reduced density matrix (A^ aa .) satisfies the following condition^: 

N 2 \ 2 max ~0(N). (36) 

Calculated wave functions and second order reduced density matrices in the vector product 
method are not non-interacting in general. Then there is a possibility that the supercon- 
ductivity of solid can be explained by the results of the vector product method. 
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